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' Based on the Relativistic Mean Field (RMF) approach the existence of 

o\ '• 

Q^ ' the broken pseudo-spin symmetry is investigated. Both spherical RMF and 

constrained deformed RMF calculations are carried out employing realistic 

s : 

Lagrangian parameters for spherical and for deformed sample nuclei. The 

m : 

quasi - degenerate pseudo-spin doublets are confirmed to exist near the fermi 

> ' 

' surface for both spherical and deformed nuclei. 
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Pseudospin symmetry has been discovered in nuclear physics nearly 30 years ago 
The recent claim jl that pseudo-spin symmetry may arise due to near equality in magni- 

. 

tude of attractive scalar and repulsive vector fields in relativistic mean theory, has revived 

m > '■ 

the activity related to the understanding of the origin of this symmetry in real nuclei. The 
concept of pseudo-spin symmetry |I|,0] is based on the experimental observation of the ex- 
istence of quasi-degenerate doublets of normal parity orbitals (n, £, j = £ + ~) and (n — 1, 
£ + 2, j — £+ |) such as (4si/ 2 , 3c/ 3 / 2 ), (3rf 5 / 2 , c 2g 7 / 2 ) etc., in the same major shell. Since for 
spherical systems the quantum numbers j n are conserved, the pseudo-spin angular momenta 
5=1/2) satisfy j =j = l±\. 

In order to interpret this near degenerate pair of j = £ + 1/2 and j = £ + 3/2 states 
as pseudo-spin doublets corresponding to rh s = ±1/2, £ has to be £ + 1. It then follows 
for the major oscillator quantum number: N = N — 1, for the radial quantum number 
n = (N — £) — 1 and for the parity n = —n. For zero pseudo-spin orbit splitting, the pseudo- 
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spin multiplet will be degenerate. Thus the pair of orbitals (4si/ 2 , 3<i 3 / 2 ) and (3<i 5 / 2 , 2g 7 / 2 ) 
can be viewed as the (2p 1 / 2 , 2p 3 / 2 ) and (I/5/2, I/7/2) pseudo-spin doublets. The symmetry 
can also be investigated in deformed nuclei. In the asymptotic Nilsson scheme one finds the 
pseudo-spin quantum numbers (N = N — 1, n 3 = n 3 , A = A + 1 and Q = fl) . Therefore, 
the Nilsson orbitals [N, n 3 ,A,Q = A + 1/2] and [N, n 3 , A + 2, Q = A + 3/2] can be viewed as 
the pseudo spin-orbit doublets [N, h 3 , A, ft — A ± 1/2] ||. 

Apart from the rather formal relabeling of quantum numbers various proposals for an 
explicit transformation from the normal scheme to the pseudo-spin scheme have been made 
in the last twenty years and several nuclear properties have been investigated in this scheme 
|6|-||. However, the origin of pseudo-spin symmetry remained unknown until the recent 
observation of Ginocchio P,|10|| where for the first time the origin of this symmetry is claimed 
to be revealed as due to the near equality in magnitude of the attractive scalar and repulsive 
vector fields in relativistic theories. Here in this letter we follow this idea and investigate 
to what extent the pseudo-spin symmetry is broken for realistic cases. For this purpose 
we concentrate as well on spherical as on deformed nuclei and we use the framework of 
relativistic mean field (RMF) theory [JTTJ] . It has been shown that this phenomenological 
approach is very successful in describing the ground state nuclear properties of spherical, 
deformed and also for nuclei far away from the beta stability line (see for example |12|-|l^1. 

The RMF starts with a Lagrangian density describing the nucleons as Dirac spinors ip, 
of mass m, interacting via the meson (er-, to-, and p-) and the electromagnetic fields. The 



standard Lagrangian density used in the RMF theory is written as 13 



C = ip (j(id - g^uj - g p pr - eA) - m - g a a) ip 
+ 1(^)2 _ U{a) _ Ifi^ + l -mlu 2 

~\^ V + \m 2 p p 2 - \F„JV (1) 
It includes a nonlinear self-interaction U{a) of the cr-field: 
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U(a) = l -m 2 y + l -g 2 a" + \g,a\ (2) 

which takes into account in a phenomenological way the density dependence of the param- 
eters of the model. m a (g a ), m^^g^), m p (g p ) are the respective meson masses (coupling 
constants) and g 2 and g% are the coupling strengths of the nonlinear sigma field U(a). 

It is straightforward to write the coupled baryon spinor and the mesons mean field 
equations. Starting from the Dirac equation for the single nucleon radial wave function with 
the spherical attractive scalar (S = —g a cr) and the repulsive vector (V = g^uS) potentials 
and following the standard procedure, by eliminating the small components (gi), the large 
components (/;) obey the following second order differential equation: 

J_V* (1 + «±1) ) f 

\ 2m- E- (S + V) \dr r )) Jl 

= _ { 2 m -E-(S + V)) (E-(S- V)) U (3) 
Here the eigenvalues denoted by n iy of the operator — /3(S • L + 1) are given by 



I) for 5i = U±\, (4) 



and 5" (V) are the derivatives of the potentials S (V ) with respect to r. The binding energy 
E > is measured with respect to the nucleon mass M in natural units h — c — 1. 

On the other hand eliminating the large component fi we have for the small component 
gi the following second order differential equation: 

J_V72 _ S '~ V ' (jL _ \ 

\ E-(S-V) \dr r j ] 9i 

= ( 2m - E - (S + V)) (E-(S- V)) 9i . (5) 

For the case of equal strengths, S = V, the Eq. (^) reduces to: 

- V 2 g t + E(S + V) 9l = E(2m - E) 9l . (6) 

Clearly Eq. @ has an energy dependent potential (E(V + S)) and has the eigenvalue 
E(2m — E). After scaling the radial variable r = x/(y/~E), the potential has a complicated 



(y/~E) dependence i.e., Sfx/VE^j + V (x/y/E). In such a situation this equation (Q) is 
no longer a normal Schrodinger eigenvalue equation. Further, it is obvious that in this 
equation all solutions with "bound" states in the Fermi sea with E > are shifted to one 
degenerate eigenvalue with E — 0, which, in fact, is not bound. The corresponding wave 
functions are not normalizable. This indeed is an unphysical situation. This equation is the 
same as the equation (3) of Ref. M in the scaled variable x when written in terms of the 
partial waves and using the relation £(£ + 1) = k,(k — 1). Here £, the angular momentum 
of the lower component gi is identified with the pseudo-spin angular momentum (£). This 
is the pseudo-spin symmetry limit of Ref. Q, where the doublets j = £ ± 1/2 with the 
same £ are degenerate. However, in this limit only the Dirac sea states exist and no Dirac 
valence bound states and therefore contradicts reality. According to these considerations in 
all realistic situations the pseudo-spin symmetry must be broken. Therefore the question 
arises, to which extent it is broken in real nuclei. So far only the spherical case has been 
investigated for square well potentials j| and for spherical solutions of the RMF equations 
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In the present letter we investigate the broken pseudo-spin symmetry both for the spheri- 
cal and deformed nuclei within the relativistic mean field approach. For our study, we choose 
208 Pb as a representative of spherical nuclei and 154 Dy as a representative of deformed nu- 
clei. We use in our calculations the Lagrangian parameter set NL3 [17]] which successfully 
reproduces the ground state properties of nuclei, spread over the entire periodic table. The 
other parameter sets like NL1 and NLSH (see |14]]) are expected to give almost identical 
results for these nuclei. 

First, spherical RMF calculations in the coordinate space are carried out for 208 Pb. The 
calculated binding energy and the charge radius agree remarkably well with the experiment. 
The calculated single particle energies for the bound orbitals near the fermi surface are shown 
in Fig. (|I]a) for neutrons and protons. It is clear from the figure that the pairs of bound 
neutron valence orbitals (2^7/2,3^5/2) and (lzii/2,2g 9 / 2 ) which correspond to pseudo-spin 
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doublets (2/7/2,2/5/2) and (l/tn/2, I/19/2) respectively, are quasi-degenerate indicating only 
a small breaking of pseudo-spin symmetry. The same is more or less true for the pairs of 
neutron hole ((2/ 5 /2,3p 3 / 2 ), (1/19/2,2/7/2)), proton valence (particle) (1/19/2,2/7/2)), and proton 
hole ((2cZ 3 /2,3si/ 2 ), (lg7/2,2cZ 5 /2)), orbitals forming the pseudo-spin doublets. But here the 
energy separation between the partners of the respective doublets is relatively larger. The 
larger is the binding energy the larger is the separation. This indicates that the concept 
of the pseudo-spin symmetry becomes better and better for the orbitals as their energies 
approach closer and closer to the continuum. This is consistent with the results found in 
Ref. [§J for the square well potentials. In addition, the energy separation becomes larger, if 
the pseudo-orbital angular momentum (£) increases. The dependence of the energy splitting 
of the pseudo-spin partners on the energy E and on the pseudo-orbital angular momentum 
£ can easily be understood from Eq.([5]). For a given pseudo-orbital angular momentum £ 
the term in Eq.(|]) which splits the pseudo-spin partners is: 

S'-V K (7) 
(S-V)-E r { ' 

It has the energy dependence (E — (S — V)) in the denominator. Now (S — V) is about 
50 MeV. Bound states in the Fermi sea have a binding energy E < 50 MeV. For increasing 
binding energy E , i.e. going to more deeply bound states, the denominator decreases. This 
then results in a larger energy splitting between the pseudo-spin partners. For example for 
the orbit £ = 3 the energy splitting between the pseudo-spin partners (1*77/2 and 2d 5 / 2 )) 
will be relatively larger as compared to that between (2^7/2 and 3^5/2)). In addition, the 
bigger is the value of I the larger is the splitting. For instance,the energy splitting between 
the pseudo-spin partners (lin/2 and 2g 9 / 2 ) corresponding to £ = 5 is relatively larger as 
compared to that between the partners (2<7 7 / 2 and 3rf 5 / 2 ) which corresponds to £ = 3, in the 
same major shell. Interestingly, the sign of the energy splittings between the partners of 
the neutron valence doublets is opposite to those of the neutron hole, proton particle and 
proton hole doublets. 

The normalized single nucleon wave functions (both large (/) and small (g) components) 
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are plotted for the pseudo-spin partners corresponding to the valence neutron pairs, the 
neutron hole pairs and valence proton pairs in Figs. [I] (b), (c) and (d) respectively. The 
phase of the lower components (g) of one of the partners is reversed while plotting, in order 
to exemplify the differences in the lower components of the pseudo-spin partners. Clearly, 
the lower components are much smaller in magnitude as expected and are almost equal in 
magnitude. In the case of exact pseudo-spin symmetry, the lower component of the pseudo- 
spin partners should be identical (except for the phase). The very small differences between 
these g's which mainly appear around the surface are negligible for the pseudo-spin partners 
having very small binding energies. 

Next we consider deformed systems and impose constraint on the quadrupole moment. 
Constrained relativistic Hartree calculations have been carried out for the nucleus 154 Dy. 
The coupled differential equations for the spinors and the meson fields are given in Ref. 
13fl . They have been solved by expanding the spinors and the meson fields in terms of 



anisotropic oscillator wave functions. Numerical details are given in Refs. |TJj and [JT8 
Pairing correlations are treated in the constant gap approximation and the Lagrangian 
parameter set NL3 |HJ is used. The calculated potential energy surface is shown in Fig. |2|. 



The value of the calculated ground state deformation parameter (3 2 is 0.202 which is to be 
compared with 0.237, the corresponding experimental value. The calculated ground state 
binding energy 1262.95 MeV differs from the corresponding experimental value by merely 
1.2 MeV. 

The energies of the bound neutron pairs of orbitals corresponding to pseudo-spin doublets 
are plotted against the deformation /3 2 ranging from 0.0 to 0.5 in Fig. || The asymptotic 
Nilsson quantum numbers [N, n^, A, Q] are good for large values of the deformation fa- The 
pseudo-spin doublets [N, n 3 ,A, Q = A ± 1/2] || are indicated by [N, n 3 ,A] f and j in 
the figure. For zero deformation (/?2 = 0) the orbitals are indicated by the corresponding 
spherical states. The figure reveals the following: 

(i) The energy splitting between the pseudo-spin partners is smaller for the valence orbitals 
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and for the partners just below the Fermi surface. 

(ii) This energy difference is relatively larger for the partners having larger pseudo-spin 
angular momentum (£). 

(iii) In general, this separation stays almost constant and does not vary with deformation 
after reasonable value of /3 2 . 

(iv) The energy difference between the j and the j partners always remains positive except 
for [404], where there is crossing at around (3 = 0.3. Such a crossing is not very 
unusual, it has also been observed in Ref. ||. 

These systematics are consistent with those observed in the spherical case above. A 
similar plot for the proton pseudo-spin doublets shown in Fig. [| reveals identical systematics 
as those observed for the neutron case (Fig. |||). It is interesting to note that in Ref. || the 
energy difference between the valence neutron pseudo-spin partners is negative (opposite 
to ours) while it has the same sign as ours for protons. This may be due to the negative 
value obtained for Vg s , the strength of the pseudo-spin orbit interaction, from the Nilsson 
parameterization for 82 < iV < 126. 

Similar calculations have also been carried out for other spherical and deformed nuclei 
and they show identical systematics. The conclusions presented here, are therefore rather 
general. 

In conclusion, it is shown in the relativistic mean field framework that quasi-degenerate 
pseudo-spin doublets do exist near the fermi surface for both spherical and deformed nuclei. 
The pseudo-spin symmetry is restored better and better as one moves closer to the continuum 
limit. These conclusions confirm the findings of Ginocchio PlJlOl . 
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FIG. 1. Pseudo-spin splitting in the spherical nucleus 208 Pb: (a) single particle spectra in the 
vicinity of the Fermi surface for neutrons v and protons (jr) and large (/) and small (g) components 
of the Dirac wave functions for the pseudo-spin doublets v2d (b), v2f (c) and ir2g (d) 

FIG. 2. Energy surface of the deformed nucleus 154 Dy as function of the quadrupole moment q 
in units of barn 

FIG. 3. Single particle energies of the deformed Dirac equation for the neutrons in the nucleus 
as a function of the quadrupole deformation parameter fa ■ Asymptotic pseudo-spin quantum 
numbers are given and the pseudo-spin partners are indicated by arrows j and [ 

FIG. 4. Single particle energies for protons in 154 Dy, for details see Fig. || 
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